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Atomic level defects such as dislocations play key roles in determining the macroscopic prop¬ 
erties of crystalline materials Their effects are important and wide-reaching, and range 
from increased chemical reactivity to enhanced mechanical properties to vastly in¬ 
creased rates of crystal growth. Dislocations have therefore been widely studied using tra¬ 
ditional techniques such as X-ray diffraction (XRD) and optical imaging. More recently, 
advances in microscopy have allowed their direct visualization. Atomic force microscopy 
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(AFM) has enabled the 2D study of single dislocations while transmission electron mi¬ 
croscopy (TEM), which was initially limited to 2D projections of thin specimens, can now 
visualize strain fields in 3D with near atomic resolution However, these techniques can¬ 
not offer in situ, 3D imaging of the formation or movement of dislocations during dynamic 
processes such as crystal growth and dissolution. Here, we describe how Bragg Coherent 
Diffraction Imaging (BCDI) can be used to visualize in 3D the entire network of dislo¬ 
cations present within an individual crystal. Using calcite (CaCOs) single crystals, we also 
use BCDI to monitor the propagation of the dislocation network during repeated growth and 
dissolution cycles, and show how this is intimately linked to the growth and dissolution mech¬ 
anisms. These investigations demonstrate the potential of BCDI for studying the mechanisms 
underlying the response of crystalline materials to external stimuli. 

Crystal growth and dissolution processes have been studied for over a eentury due to their 
signifieanee to fields sueh as geology, eorrosion, eatalysis and the synthesis of nano-struetures, 
through a desire to understand the link between mieroseopie and maeroseopie proeesses, and of 
eourse due to our innate faseination with sueh beautiful struetures. These investigations eonfirmed 
that dissolution and growth proeeed by analogous meehanisms and also identified that crystal¬ 
lographic defects, and in particular screw disloeations, are of fundamental importance to crystal 
growth and dissolution processes!^. This can be attributed to the faet that screw disloeations eause 
deformation (and therefore strain) of the adjaeent erystal lattiee, whieh in turn alters the aetivation 
barrier for growth and dissolution in the vieinity of the disloeation A full pieture of erystal 
growth and dissolution meehanisms ean therefore only be obtained by studying the relationship 
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between the evolution of the network of dislocations within a crystal and its morphology during 
these processes. 

In this article, we demonstrate how BCDI can be used to study the role of dislocations in 
dictating the mechanism of growth and dissolution of calcite crystals. Calcite was selected for 
study as it is one of the most-studied inorganic crystals, and the crystals readily grow to a few 
microns in size, with well-defined morphologies. BCDI is an imaging technique that uses coherent 
X-rays to image the density (and morphology) of a crystal, and importantly, the strain within it 
fHESl Illumination of a crystal that is smaller than the coherence volume of the beam generates a 
coherent X-ray diffraction (CXD) pattern due to scattering from all parts of the crystal. An image 
of the crystal morphology can then be generated from the CXD pattern, where the phase of the scat¬ 
tered wave is reconstructed using iterative phase retrieval algorithms EnnEHSI reconstructed 
density is complex valued, with the amplitude containing information about the electron density, 
p(r). Phase shifts in the reconstructed complex density arise from strain (internal deformation) in 
the crystal lattice. The phase is proportional to the vector displacement field, u{r) of the atoms 
from the ideal lattice points and the scattering vector Q via 4>{r) = Q ■ u{r) (see supplemen¬ 
tary information for further details). For a single Bragg peak a single projection of u is obtained 
(projected displacement) and components of u perpendicular to Q will not be observed. It is this 
sensitivity to deformations that makes BCDI ideal for studying defects within crystals. 

Calcite crystals were precipitated by placing 100 /iL droplets of a solution containing CaCl 2 , urea 
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and urease on hydroxyl terminated, self-assembled monolayers (SAMs) supported on gold thin 
films. CaCOs preeipitation then oeeurred on enzymatic hydrolysis of the urea to ammonium and 
carbonate This method was selected as it generated a high density of {104} oriented calcite 
rhombohedra with average diameters of ~ 1.25 /xml^. BCDI experiments were then carried out at 
beamline 116 at Diamond Light Source. Individual calcite crystals below 2 /xm in size were illumi¬ 
nated with monochromatic, 8 keV X-rays, and diffraction was recorded at the {104} Bragg peak. 
3D diffraction data sets (Supplementary Information Fig. 1) were obtained by rocking an isolated 
calcite crystal through its Bragg peak, and the same crystal was monitored while undergoing cycles 
of growth and dissolution. Dissolution was achieved by depositing dilute acetic acid solution on the 
crystal, while growth was achieved by adding a drop of calcium bicarbonate solution. Alignment 
of the crystal was maintained throughout, as the X-rays were nominally unfocussed and defined 
by slits with a square opening of 200 /rm placed 0.3 m before the sample. This large beam size 
relative to crystal size also ensured that the sample was coherently illuminated 1^. 

Figure 1 shows 3D images of the crystals as iso-surface renderings of the reconstructed ampli¬ 
tudes (electron density, see supplementary information for details) (a) and phase (projected lattice 
displacement) (b) of the initial crystal (i), after growth (ii) and after two consecutive dissolution 
steps (iii and iv). The initial crystal (i) is the expected rhombohedron, where this is consistent with 
scanning electron micrographs (Supplementary Information Fig. 3). Growth of this crystal (i to ii) 
leads to an increase in size and a smoothing of the faces exposed to the bulk solution, while the 
face in contact with the SAM remains unchanged. Interestingly, two of the faces directed into the 
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solution (indicated with blue arrows) grow more rapidly than the other three. This is immediately 
indicative of a non-uniform distribution of defects, as all faces of an entirely perfect rhombohedron 
would be expected to grow at the same rate. Images of the corresponding projected displacements 
are shown in Figure lb, where this is mapped onto an iso-surface with red and blue representing 
lattice contraction or expansion respectively by half a lattice spacing. Comparison of the projected 
displacements before (i) and after (ii) crystal growth shows that these do not grow significantly 
with the crystal but remain maximal at the edges. This is indicative of the presence of active 
growth fronts 1^. 

The dissolution steps (Fig. 1, ii to iii and iii to iv) show that the crystal faces retreat along all 
directions, but that this is more pronounced at certain sites (indicated by red arrows). Initial signs 
of changes in the crystal shape and the onset of etch-pit formation are visible after the first stage 
of dissolution (iii), which leads to an increase in the specific surface area and roughness of the 
crystal. The etch pits are also associated with higher levels of deformation/strain. That relatively 
little change occurs in the crystal face adjacent to the substrate is consistent with the intimate con¬ 
tact of the SAM with this nucleation face. The second dissolution step (Fig. 1 (iii-iv) and Fig. 5 
(iii-iv)) results in a significant change in the crystal morphology, and the production of a porous 
isometric form ^ that is quite distinct from the original shape and which can be attributed to the 
removal of defect outcrops at the crystal surface and the coinciding etch pits formed 1^^^. Looking 
in turn at the lattice deformation, it is evident that strain present at the crystal surface reduces with 
the increased dissolution. This indicates that the least stable (more strained) regions dissolve first. 
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leaving behind a more stable eore (Supplementary Movies 1-4). 


The projeeted displaeement images also reveal a further intriguing feature, whieh is indieated 
by the grey arrows in Fig. lb, and whieh is present throughout the growth and dissolution of the 
erystal. This region possesses both a hollow eore and a spiral phase, where this eombination of 
features is eharaeteristie of serew disloeations; a number are highlighted in Fig. 2. Disloeations are 
eharaeterized by a Burgers veetor that measures the topologieal shift of the erystal along the dislo- 
eation line, where this is usually a lattiee veetor of the erystal 1^. Whenever there is a eomponent 
of the Burgers veetor parallel to the disloeation line, it has a serew disloeation eharaeter eausing 
the lattiee to spiral around the disloeation; in this way erystal growth (and dissolution) ean be fa- 
eilitated. Confirmation that this feature indeed eorresponds to a serew disloeation was obtained by 
reeording the polar angle dependenee of the displaeement assoeiated with the eore (indieated by 
the eirele in Figure 2a) over the growth/ dissolution eyele of the erystal (Figure 3a). An approxi¬ 
mately linear relationship was observed and is eonsistent with what is predieted by linear elastieity 
theory (see supplementary information). To further eonfirm the nature of the identified disloea¬ 
tion, a eomparison is provided between a simulated serew disloeation and simulated serew after 
BCDI proeessing. The model serew disloeation used in the simulation is shown at atomie resolu¬ 
tion in Figure 4a, with the resulting displaeement of atoms from their ideal lattiee positions given 
in Figure 4e. The BCDI experiment proeessed simulated serew disloeation is shown in Figures 
4b and 4d (see supplementary information for details). Figures 4e and 4f provide a eomparison of 
BCDI reeonstrueted displaeement and simulated displaeement, viewed along the disloeation line. 
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and clearly show that the low density core and spiral displacement are well-preserved after BCDI 
processing. 

Further examination of regions that showed a spiral deformation and low-amplitude eore enabled 
many additional dislocations to be identified within the imaged erystal (see supplementary infor¬ 
mation). These have been rendered and are shown in Figure 5 and Supplementary Movies 5-8. 
The initial crystal (i) possesses several dislocations which are located relatively close to the erystal 
surfaee (down to 200 nm) and are found predominantly parallel to {104} planes. These disloea- 
tions are mainly assoeiated with the faees exposed to the bulk solution, and are visible throughout 
a single crystal growth and dissolution cycle (steps i-iii), during whieh they (ii) increase and (iii) 
deerease in length. The reprodueibility of locating the same dislocation across independently re- 
eonstrueted data sets provides further eredenee to the reeovered images. 

Disloeations are often found to oeeur in loops or pairs with opposite Burgers veetor, as this re- 
duees the lattiee potential energy and resulting long-range strain assoeiated with these features. 
They also are stabilized near crystal surfaces, as found in this work, such that dislocation motion 
can transport material into and out of erystals from the growth solutionis. As intuitively expeeted, 
the most rapidly growing erystal faees were observed to have the highest number of surfaee dislo¬ 
eations. The distribution of these defeets within the erystal therefore plays a part in determining 
the morphology of the produet erystal. Considering the dissolution process, the loeations of the 
disloeations eoineide with the position of the eteh pits that appear at the crystal surface during 
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dissolution. This suggests that the developing pits follow the cores of the dislocations, where this 
effect can be attributed to a reduced activation barrier to dissolution due to stored strain energy 
in the defects 14. New dislocations could be identified during the dissolution (iii to iv), and by 
the final stage (iv) many of the original dislocations had been annihilated. The loss of faceting 
and dislocations near the surface support the hypothesis that the least stable regions dissolve first, 
leaving behind a more stable core. 

To determine the overall effect that the growth/ dissolution had on the strain and deformation 
of the crystal, the gradient of the displacement was calculated and its magnitude was plotted as a 
function of fractional size (Figure 3b and c). The magnitude of the gradient was calculated over 
successively larger shells and averaged over all directions and is plotted as a % of deformation 
(strain) relative to the lattice constant for the 104 reflection (see supplementary information for de¬ 
tails) for the directions parallel and perpendicular to the scattering vector (which is approximately 
parallel with y). What is evident (Fig. 3b and c) is that the initial crystallite experiences an increas¬ 
ing strain from the center, outwards. This can be compared to a crystal formed after the first growth 
stage (ii), which shows an overall flatter initial strain, and then increases more rapidly above a frac¬ 
tional size ~ 0.6. The smoother transition and lower strain for the direction parallel suggest some 
relaxation of the crystal after the initial growth and can also be seen in the slices in Fig. 2. These 
data are therefore consistent with the assertion that the effect of the surface only penetrates to a 
finite distance into a crystal. With the increased roughness and etch-pit formation that occurs in 
the first dissolution step (iii), the overall strain is higher, demonstrating that these surface features 
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affect the entire crystal, leading to an increase in overall deformation. After the final dissolution 
(iv) when the majority of surface dislocations have been removed, the strain becomes almost flat. 
This indicates that surface effects are less pronounced, despite the increased relative surface area. 

With its ability to simultaneously generate 3D images of the strain within a crystal and the gross 
crystal morphology, BCDI provides an extremely powerful way of visualizing the network of dis¬ 
locations present within an individual crystal. Importantly, BCDI can be performed in the absence 
of any sample preparation, which means that it can be used to elucidate the effects of dislocation 
networks on dynamic crystal behavior. Using the example of crystal growth and dissolution, we 
have located dislocations within calcite single crystals, and then demonstrated that their distri¬ 
bution within the crystal dictates rapidly-growing directions. Conversely, preferential dissolution 
and etch pit formation was also observed within the vicinity of the dislocations, which provides 
evidence that the stored energy within the dislocations affects the crystal response. The ability to 
view these fundamental processes using BCDI represents an important step forward in elucidating 
the nanoscale mechanisms underlying crystallization processes. Importantly, BCDI also opens the 
door to 3D visualization of the role of dislocations in the response of a crystal to a huge range of 
experimental conditions, such as temperature changes or mechanical force. 
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Figure 1 : Growth and dissolution of calcite observed by BCDI Shown is an iso-surface render¬ 
ing of the (a) electron density (reconstructed amplitude) and (b) projected displacement (phase) 
from (i) initially deposited calcite crystal, (ii) after secondary growth, and (iii & iv) after consec¬ 
utive dissolution steps. The scale bar is 1 /xm. Three different viewing angles of the crystal are 
shown in (a) and (b). Given are top down (left), side (middle panel) and bottom up (right) perspec¬ 
tives of the imaged calcite crystal sitting flat on the substrate surface. Perspectives highlight the 
shape transition that occurs during growth (i-ii) and dissolution (ii-iv). Prominent surface advance 
(growth, blue arrows) and retreat (dissolution, red arrows) directions are shown. The gray arrow 
points towards the primary dislocation continuously identifiable during crystal growth and disso¬ 
lution (i-iii). The beam direction is along the z axis, with the y axis vertical. The sample/substrate 
is located at a set scattering angle towards the beam direction (z) and spanned plane (x-z). 
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Figure 2: Slices showing the electron density and projected displacement during growth and 
dissolution Shown is the initial crystal (i), after growth (ii) and repetitive dissolution steps (iii & 
iv) for two viewing directions (a) and (b). The amplitude is shown in the top row (black back¬ 
ground) with the dislocations highlighted by white arrows featuring a low-amplitude core. The 
phase presented in the bottom row also shows the present, selected dislocations highlighted by the 
dark blue arrows. Particularly evident is the spiral phase (displacement) that is characteristic of a 
screw dislocation. The iso-surface to the right shows the location of the cut-planes. The scale bar 
is 1 /im. 
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Figure 3: Displacement and strain line plots. Plotted in (a) is the recovered projected displace¬ 
ment measured at development stages (i-iv) as a function of theta for the dark blue circles with 
arrows from Fig. 2A. This is compared to the displacement expected by linear elasticity theory 
(LET), (b) The strain component parallel and perpendicular (c) to the scattering vector over the 
growth/dissolution cycle (i-iv) plotted as a function of fractional crystal size (center of the crystal 
0 to crystal surfaces 1). This graph highlights the diminishing relevance of surface effects with (ii) 
growth of the crystal, and its increasing relevance with (iii & iv) dissolution. The strain is calcu¬ 
lated as the magnitude of the gradient of the displacement and averaged over successively larger 
shells. 
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Figure 4: Simulation of a screw dislocation Iso-surface rendering of a screw dislocation (a) with 
atomic resolution, and by filtering (b) the Fourier transform of (a) with a Gaussian centred on 
the (0,1,0) Bragg peak, replicating the BCDI experiment. The displacement is rendered onto the 
iso-surface for the atomic resolution (c) and the phase is rendered onto the surface for the BCDI 
simulation (d) showing the spiral phase centred around the dislocation core, (e) and (f) show 
another view of the dislocation for the atomic resolution (e) and BCDI experiment (f), revealing 
the low density region at the core in the BCDI experiment simulation. It should be noted that for the 
atomic resolution the displacement is mapped onto the iso-surface where for the BCDI example, 
the resultant phase is mapped on to the iso-surface. 
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Figure 5: Iso-surface rendering of defect network within a calcite crystal. Iso-surface ren¬ 
derings of dislocations present within the crystal imaged (i) before, and after (ii) growth / (iii-iv) 
dissolution. The scale bar is 1 /rm. The evolution of the dislocations is evident through crystal 
growth (i to ii) and dissolution (ii-iii, iii-iv), where these are shown from left to right. Three dif¬ 
ferent viewing angles of the crystals are provided, (a) from base, (b) from top and (c) from side. 
Dislocations are predominantly located near the surface and on the faces that grow most rapidly. 
The phase has been mapped to the iso-surfaces of present dislocations, showing the characteristic 
spiral. 
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Methods 


Materials and General Preparative Methods Analytical grade, CaCl2-2H20, Acetic Acid, Urea 
and Urease (canavalia ensiformis, subunit molecular weight: ~90.770 kDa) were purchased from 
Sigma-Aldrich and were used as received. Aqueous solutions were prepared using Milli-Q Stan¬ 
dard 18.2 MUcm. Solutions of 11-Mercapto-l-undecanol (Sigma-Aldrich) were freshly prepared 
in laboratory grade ethanol, and experiments were performed at a constant temperature of 21°C. 
Reagent containing glassware was soaked overnight in 10% w/v NaOH, followed by rinsing with 
dilute HCl and washing with Milli-Q water. 

Substrate Preparation Glass slides and crystallizing dishes were placed overnight in Piranha 
solution (70:30 %wt. sulfuric acid: hydrogen peroxide) and were then washed copiously with 
Milli-Q water before being exposed to the crystallisation solution. Functionalized self-assembled 
monolayers (SAMs) were prepared on freshly deposited noble metal films. Thin films were de¬ 
posited either on silicon wafers using a Mantis Qprep 250 deposition system or a Cressington 308R 
coating system at a base pressure below 10“® mbar. 2 nm of Cr were initially deposited to promote 
substrate adhesion, followed by the evaporation of 50 nm of Au (Goodfellow, 99.99%) at <0.1 
nm/s. SAMs were then prepared by immersing a prepared metal substrate in 1 mM thiol solution 
in ethanol, at room temperature for 24 h in the dark. The SAMS were then thoroughly rinsed with 
ethanol and Milli-Q water and were dried under nitrogen before usage. 

Mineral Deposition Preferentially oriented calcite was obtained by adding 350 jA urease (1 mg 
mL“^) to 1 mL of 5 mM CaCl 2 / 20 mM urea. Separate droplets (100 ^uL) of the resulting aqueous 


19 



solution were subsequently plaeed on the prepared SAM, and stored at 100% r.h.. Hydroxyl- 
terminated SAMs were used as the erystallization substrate to obtain ealeite rhombohedra that 
were prineipally {104} oriented. Enzymatie hydrolysis of urea to ammonium and earbonate ereates 
the required supersaturation profile, where the gradual inerease in supersaturation ensures that a 
suffieient density of single erystal partieles of CaCOa is obtained. 

Overgrowth of Initial Deposits Observations of erystal growth were earried out on samples 
whieh had previously been eentered into the Bragg eondition. A 50 /i/ volume of ~ 1 mM ealeium 
biearbonate solution was pipetted onto the sample, where evaporation and CO 2 out-gassing results 
in a supersaturation inerease whieh eauses erystal growth. After evaporation of the droplet, the 
seleeted erystals were then re-analyzed. Caleium biearbonate solutions were prepared aeeording 
to the method of Kitano 1^. 100 mg of CaCOa were added to one litre of Milli-Q water. C02{c,) was 
then bubbled through this solution for three hours followed by filtering through a 200 nm Isopore 
GTTP membrane filter (Millipore) resulting in an initial pH of 6.4. 

Partial Dissolution Dissolution of BCDI imaged ealeite erystals was aehieved by depositing 50 
/xL of 0.1 wt% aeetie aeid solution onto a 0.4 em^ substrate with an estimated number density of 
0.1 erystals/ /rm^ with an average size of 1 /im. The solution was then removed after 60 seeonds. 
This was followed by addition and removal of a drop of ethanol to wash the sample, and a further 
diffraetion pattern was eolleeted. The proeess was then repeated between two and three times to 
obtain diffraetion data for sueeessive dissolution stages of the initially imaged single erystal. 
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Enzymatic Hydrolysis (Urea-Urease) catalyzed Precipitation Calcium carbonate deposition 
using the enzymatic hydrolysis of urea by urease, and droplets placed on the substrate, allows 
controlled precipitation of crystals independent of droplet wetting behavior and gas-liquid inter¬ 
face area By variation of the urea and urease concentrations, which determines the rate of 
production of carbonate and ammonium, the rate of production of CaCOs can be controlled. A 
slow increase in supersaturation was required to give crystals that were principally unstrained and 
which were present in a high enough number density to enable crystals to be easily located on the 
substrate during BCDI imaging. The key underlying reactions are presented below. 


{{NH2)2C0) + 2 H 2 O 2{NH+) + {COI-) 

(1) 

{{NH 2 ) 2 C 0 ) + H 2 O ^ 2(iVff3) + {CO 2 ) 

(2) 

{NH 3 ) + H 2 O ^ {NHt) + OH- 

(3) 

{CO 2 ) + H 2 O 0 {H 2 C 0 ^) 

(4) 

{H2C0^) ^ {HCO^) + {H+) 

(5) 

{HCO^) 0 {COI-) + (iT+) 

(6) 

{COI-) + {Ca^+) 0 CaCOs 

(7) 


External Characterization CaC03 precipitates at different growth and dissolution stages were 
characterised by Scanning Electron Microscopy (SEM). Electron micrographs of uncoated spec¬ 
imen were obtained using a EEI Nova NanoSEM 650. Crystal growth and dissolution rates of 
outgrown individual calcite crystal were obtained using an inverted Olympus IX-70 confocal mi- 
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croscope. 


Experimental Setup BCDI experiments were performed at Beamline 1-16 at the Diamond Light 
Souree (DLS) UK. An undulator produeed X-rays whieh were monoehromatized to an energy 
of 8 keV using a ehannel-eut Si(lll) monoehromator. Caleite erystals supported on substrates 
were plaeed on a diffraetometer whieh had its rotation eentre aligned with the X-ray beam. Slits 
were used to reduee the area illuminated by the X-rays. An X-ray sensitive, CMOS pixel detee- 
tor (MedipixS ) with 256x256 square pixels of side length 55 //m was positioned at the desired 
diffraetion angle for an off-speeular {104} refleetion at a distanee of 2.7 m from the sample. An 
evaeuated flight tube was plaeed between the deteetor and sample to reduee air absorption and seat- 
ter. To measure its full 3D diffraetion pattern, the erystal was rotated by 0.3 degree with a 0.003 
degree step size. A two-dimensional sliee of the 3D far-field diffraetion pattern was reeorded at 
eaeh angle of rotation. By staeking all of these two-dimensional diffraetion frames together, a 
eomplete 3D diffraetion pattern was obtained, from whieh real-spaee images ean be reeonstrueted 
(see Phase Retreival). Due to the small size of the erystals, ~ 1-2 //m, the illumination ean be eon- 
sidered to be almost eompletely eoherent. To ensure that the same erystal eould be analysed over 
a growth/ dissolution eyele, a loeal seareh in reeiproeal spaee was performed to piek an isolated 
single erystal. This ensured that there was little ehanee of a signal originating from another erystal 
overlapping the original erystal diffraetion. 

Data preparation For eaeh growth/dissolution stage, the erystal was rotated through its rooking 
ourve where a series of 2D diffraetion patterns were eolleoted. A total of 121 patterns were eol- 
leoted along the rooking ourve with a separation of 0.003 degrees. The exposure time for eaeh 


22 



pattern in the rocking curve was 1 second. This process (collecting 2D diffraction along the rock¬ 
ing curve) was repeated three times for each growth/dissolution stage. This resulted in three data 
sets for each stage which were then summed together to form a single 3D diffraction pattern which 
can be seen in Supplementary Information Fig. 1. The alignment of the three sets was done by 
calculating their centre of masses and shifting the data sets to the centre of mass of the first data set 
(to the nearest pixel). This data set was then thresholded, where all pixels with 2 or less photons 
were set to 0. Several pixels on the detector were reading high values even when no photons were 
present (‘hot pixels’). These pixels had their values set to 0. In each of the two detector pixel direc¬ 
tions, the data were binned by a factor of 2 and cropped to a final size of 80 pixels while in the scan 
direction the data was not binned and was padded (with zeros) to a final size of 128 pixels. After 
the initial growth stage, diffraction from another (different) crystal was observed on the detector. 
This diffraction was identified as coming from a different crystallite by the fact it was spatially sep¬ 
arated on the detector (from the original diffraction) with no visible interference fringes between 
this diffraction and the original. This diffraction was removed by setting it to 0. This summed, 
thresholded, zeroed and binned/padded data set was then used for the reconstruction. 
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